Revisiting the production of charmonium plus a light meson at PANDA 
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In this work, we calculate the total cross sections and the center-of-mass frame angular distributions of the 
charmonium production plus a light meson by the low energy pp interaction. The results of pp — > n'^^ with and 
without form factor (FF) indicate that the FF contribution in the calculation cannot be ignored. The obtained 
cross section of pp — > rf'JjiJ/ with FF can fit the E760 data well. We also predict the total cross sections and 
the center-of-mass frame angular distributions of pp — > oj^, which show that these physical quantities are 
dependent on Pauli (g^) and Dirac (/c^) coupling constants of the ppu) interaction. Thus, pp — > w'P can be as 
the ideal channel to test the different theoretical values of g^ and k^^. Applying the formulae of pp — > ifi'V and 
pp — > w*?, we predict the total cross sections of the pp — » rj^ and pp — » pT reactions. Our results show a 
common behavior of the chaiTnonium production with a light meson by the pp interaction, where the total cross 
section of the rj^ production is the largest one among all discussed processes. The above observations can be 
directly tested at the forthcoming PANDA experiment. 

PACS numbers: 14.40.Pq, 13.75.Cs, 13.75.Gx 



I. INTRODUCTION 

Since the discovery of Jj^ in 1974 [1, 2], more and more 
charmonia have been reported by experiment [3], which pro- 
vide an ideal platform to improve our understanding of non- 
perturbative Quantum Chromodynamics (QCD) dynamics. In 
the past decade, the observations of a series of charmonium- 
like states named as XYZ state have stimulated extensive inter- 
est in studying higher charmonia for both theorists and exper- 
imentalists. As one of the forthcoming experiments relevant 
to the study of hadron physics, Antiproton Annihilations at 
Darmstadt (PANDA) experiment at the Facility for Antipro- 
ton and Ion Research (FAIR) can serve as the investigation 
of charmonium, which is also one of the main physical aims 
of PANDA [4]. Thus, theoretically studying the charmonium 
production by the low energy pp interaction becomes an im- 
portant and interesting research work, which can give valuable 
suggestions for the forthcoming PANDA experiment. 

In Ref. [5], Gaillard and Maiani firstly calculated the differ- 
ential cross section of the charmonium production accompa- 
nied by a soft pion in the low energy pp interaction, where two 
hadron-level diagrams were introduced by the Born approxi- 
mation. They indicated that the corresponding cross section is 
proportional to the partial decay width of charmonium decay 
into pp [5]. The authors of Ref. [6] further studied the cross 
sections of the charmonium (^) production plus a light meson 
(m) by the pp ^ 'I' + m processes, which can be related to the 
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measured partial decay width of charmonium decay into ppm 
[6]. By this approach, the cross sections of pp — > J/ifr{t//')m 
(m = 7r°, ?7,p°, 0), Tj', (p) are estimated. Among these predicted 
cross sections, the cross section of the pp — > jf'JIkjj process 
can reach up to 300 pb when taking yfs - 3 GeV, which is 
close to the E760 data taking V^ = 3.5 ~ 3.6 GeV [7]. Later, 
Barnes and Li developed the initial state light meson emission 
model, which was applied to study the near threshold associ- 
ated charmonium production process pp — > n'^^', where 'I' de- 
notes rjc, Jjip, ^',XcO,Xci [8]. By the initial state light meson 
emission model, they calculated the differential cross section 
and total cross sections of pp — > n'^^' processes. In addition, 
{dcr/dQ.}, the center-of-mass frame unpolarized angular dis- 
tribution, was predicted for pp — > ;r"*I' [8]. In Ref. [9], Barnes 
et al. further indicated that the cross section of pp — > jfJIij/ 
near threshold may be aff'ected by the Pauli Jjif/pp coupling, 
which will be an interesting research topic in PANDA [9]. 

By the initial state light meson emission model [8], pp — > 
tt"*!' occurs via a proton exchange between p and p. In 
Ref. [8], authors treated the ppm and Jl^ipp couplings as the 
point-like interaction vertices. However, in reality we should 
consider the structure effect of the ppm and J/t/^pp interac- 
tions. Thus, for reflecting such structure effect, the form fac- 
tor should be introduced in the ppm and J/ifrpp interaction 
vertices, which was listed as one of the future developments 
of the initial state light meson emission model [8]. Along this 
way, in this work we revisit the production of charmonium 
plus a light meson in the low energy pp interaction by consid- 
ering the contribution of form factor (FF) to these processes. 
The comparison of the results with and without including FF 
in the calculation can reveal the difference under two cases, 
which will be tested at the forthcoming PANDA experiment. 
By this study, we can not only extract the physical picture de- 



picting the production of charmonium plus a light meson in 
the low energy pp, but also learn what kind of FF to be suit- 
able to describe the structure effect of the ppm and J/iffpp 
interactions. These processes discussed in this work also in- 
clude the charmonium production with a light vector meson 
(a>, p), where the ppco or ppp interaction relates to both Dirac 
and Pauli couplings. At present, the coupling constants of 
Dirac and Pauli couplings of ppoj and ppp are determined by 
some theoretical groups by different processes and different 
models [10-18]. Thus, we adopt these determined coupling 
constants in our calculation, which includes the total cross 
sections, the differential cross sections {dcr/dH}. These stud- 
ies can serve as further experimental test of these coupling 
constants by the charmonium production with a light vector 
meson at PANDA. 

This work is organized as follows. After introduction, we 
present the calculation of the production of charmonium plus 
a light meson in the low energy pp interaction. In Sec. Ill, the 
numerical results are given. The last section is the discussion 
and conclusion. 



II. THE PRODUCTION OF CHARMONIUM 

As depicted by Fig. 1 , the charmonium CP) production plus 
a light meson (m) by the low energy pp interaction can occur 
via the transition of pp into T + m by exchanging a proton. 
Thus, there exist two hadron-level diagrams [5, 8] shown in 
Fig. 1 if only considering the tree-level contributions. 




FIG. 1: The diagrams describing the production charmonium plus a 
light meson by the pp interaction. Here, "P denotes charmonium (rjc, 
J/if/, if/', Xco, Xci) while m is light meson in'\ rj, p, w). 



When deducing the corresponding production amplitude, 
we use effective Lagrangian approach. The interaction ver- 
tices of ppm and pp^' include 
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where (p denotes the field of nucleon and the t is Pauli ma- 
trix. (T^iv - i[y^„'Yv]/2 and Ff^y - df^ajy - dyOjft. For the ppcj 



interaction, there exists two independent coupling constants, 
i.e, go, = gNN(d and k^ corresponding to Dirac and Pauli terms 
respectively. 

In the following, we illustrate the calculation of the char- 
monium production. Without introducing FF in each of inter- 
action vertices in Fig. 1, the resulting amplitude of pp — > n'-^^' 
is [8] 
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where nip is the mass of proton, pi, pn and k are the four 
momenta of proton, antiproton and the emitted light meson. 
s, t and M are the Mandelstam variables with the definitions 
s = O^i + pif, t^ipi- kf and u^(k- pif. g„ = gpp^ and 
g^ = gpp'f denote the coupling constants of the ppn and pp^ 
interactions, respectively, vp and Up denote the spinors of an- 
tiproton and proton, respectively. In Eq. (1), Ti is defined as 
75, -;, -/y^e^ and -ijpjs^^,, corresponding to pp -^ jt"*!* 
processes with 'I' taken as rj^, xm, Jli'i^') ^^^ Xc\^ respec- 
tively. 

For the pp — > w*P process discussed here, its production 
amplitude can be written as 
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where e* is the polarization four-vector of the emitted <jj me- 
son. The definition of T2 is the same as that of Y\ in calculat- 
ing pp — > 7r°*I'. 

Just discussed above, the FF contribution to pp — > tt"*!' and 
pp — » (jj^ is not included in Eqs. (1) and (2). Since the ppm 
and pp^ vertices are not the point-like interactions, we need 
to introduce FF in each of interaction vertices, where the FF 
not only reflects the structure effect of the interaction vertex 
but also plays an important role to compensate the off-shell ef- 
fect of the exchanged proton. For comparing the results with 
and without FF, we calculate the pp — > ifi^ and pp — » uV 
processes considering the FF contribution, where the corre- 



spending amplitudes can be expressed as 



K^^^ 



gng^Vpdn, S2) 



it - mj,) 



{u - m-p) 



Upipusi) (3) 



and 






gojg^VpiPl, So) 

(f-pT+ nip) 

X -^^-^ 

(m - mj,) 

XVpiPl, S2) 



i, - It + m,,) T T 
it - mj) 



^iT'iqt) 



Up(pusi)e* +i 



f^ajgcjg^ 



■ m,,) 



(t - mp 



(^^kyT^qf) + (T^'ky 



^iJtzlil^r^TM;) 

(u - m-) 



Up(pi,si)e* 



(4) 



where 'F^(qf) denotes the introduced monopole FF with the 
definition 'Fiqf) - (A^ - inf)/{A~ - qj). m, and qi are the 
mass and the four-momentum of the exchanged proton, re- 
spectively. A is a free parameter, which is expected to be 
around 1 GeV. In the next section, we will discuss the value 
of A adopted in our calculation. In the above expressions, q, 
and qu denote the four-momenta of the exchanged protons of 
t-channel and u-channel of the pp — > m^ process just shown 
in Fig. 1, respectively. In addition, the superscript FF is in- 
troduced for distinguishing the amplitudes with and without 
FF. As indicated in Ref. [20], the pp — > m^ processes may 
include transition distribution amplitudes [19], which gener- 
alize the form factors that we include in the ppm vertex [20] . 

We notice that the ampUtudes listed in Eqs. (1) and (2) are 
indeed transverse. When introducing FF in these amplitudes, 
we cannot make these amplitudes keep transverse. Of course, 
finding a more suitable form of FF is an interesting research 
topic, where this FF not only reflects the realistic physical pic- 
ture but also can make the corresponding amplitudes still be 
transverse. In this work, we still choose the monopole FF in 
our calculation for reflecting the realistic physical picture, and 
estimate the production of charmonium by the pp interaction 
processes. 

The general differential cross section of pp — > m^ is given 
by [3] 
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where pi„„ is the three-momentum of proton in the center-of- 
mass frame of pp. The overline indicates the average over the 
polarizations of pjp in the initial state and the sum over the 
polarization of ml'V in the final state. 

With these obtained amplitudes listed in Eqs. (l)-(4), we 
finally get the expressions of the differential cross sections of 
the pp — > w*!* processes without including the FF contribu- 
tion, which are collected in appendix A. The corresponding 



total cross sections of pp -^ w*? are shown in appendix B. 
We also confirm the deductions of the differential and total 
cross sections of pp — » n^^ in Ref. [8]. Since the formulae 
of the differential and total cross sections of pp — > n^^ and 
pp — > (jjV with the FF contribution are very complicated, we 
do not show their concrete expressions in detail, but directly 
apply their formulae to the numerical calculation. 



III. NUMERICAL RESULTS 

Before presenting the numerical result, we first introduce 
the coupling constants adopted in our calculation. For the 
processes pp — > n^^, the coupling constant of the ppn in- 
teraction was given in many theoretical works, where we 
take gNNn - 13.5 [21]. The coupling constants of char- 
monium interacting with nucleons are not well established 
at present. Thus, in this work we adopted the same val- 
ues as those in Ref. [8], where these coupling constants are 
gpp,,^ = (19.0 + 3.2) X 10-3, gppji^ = (1.62 + 0.03) x 10-^, 
gpp^, = (0.97 ± 0.04) X 10-3, gpp^^^ = (5.42 + 0.37) x 10-^, 
gppxcx - (1-03 ± 0.07) X 10 3, which were estimated by the 
measured partial width of *P ^ /?/? [3]. 

For the discussed processes pp — > w*?, there are two strong 
coupling constants g^ and k^ for the ppoj vertex. At present, 
different models and different approaches gave various the- 
oretical values for these two coupling constants, which are 
listed in Table. I. 

TABLE I: The estimated values of the coupling constants g„ and k^. 



Mechanism/Model 


Sijj 


K,o 


Paris [12, 13] 


n.i 


-0.12 


Nijmegen [14] 


12.5 


+0.66 


Bonn [15] 


15.9 





Pion photoproduction [16] 


7- 10.5 





Nucleon EM form factors [17] 


20.86 + 0.25 


-0.16 ±0.01 


QCD sum rule [18] 


18 ±8 


0.8 ± 0.4 


^Po quark model [11] 


- 


-3/2 




23 ±3 





Light meson emission model [10] 








14.6 ±2.0 


-3/2 



In addition, the masses of the hadrons involved in our 
calculation include m,fi - 135.0 MeV, m^ - 782.7 MeV, 
nip = 938.3 MeV, m,,^. = 2980.3 MeV, my/^ = 3096.9 MeV, 
m^^^ = 3414.8 MeV, m^^, = 3510.7 MeV and m^, = 3686.1 
MeV, which are from Particle Data Group [3]. 



A. The total cross sections and angular distributions of 



pp- 



Tl^^V 



Firstly, in this work we confirm the results of the total cross 
sections of pp — > tt"*!' in Ref. [8], where the pp^ and ppn 



vertices are as point-like structures. Just shown in the left- 
hand diagram of Fig. 2, the dependence of the total cross sec- 
tions of pp -^ jt'^^' on the center-of-mass energy £„„ is pre- 
sented with a comparison with the experimental data from the 
E760 experiment [7], where the cross section of pp — > jiPJIxj/ 
near 3.5 GeV was measured at Fermilab. The theoretical re- 
sult of pp — » jfij/if/ taking £„„ = 3.5 GeV is larger than the 
experimental value obviously. 
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FIG. 2: (color online). The obtained total cross section of pp — » n^^ 
and the comparison of the experimental and theoretical results. Here, 
the red points with en'ors are the experimental measurement from 
E760 [7]. The left-hand and right-hand diagrams correspond to the 
theoretical results without and with the FF contribution to pp — > 
tP^, respectively. 



In Sec. II, we investigate the structure effect of the pp^ 
and ppuj vertices on the total cross sections of pp -^ n^'''^, 
which is presented in the right-hand diagram in Fig. 2. When 
taking A = 1.9 GeV, we obtain the total cross section of 
pp -^ n^^Jlijj at Ecm = 3.5 GeV consistent with the experi- 
mental data, which indicates that the FF contribution cannot 
be ignored in studying pp — » ttP^'. Adopting the same A 
value, we also obtain the total cross sections of other charmo- 
nium productions, which are listed in the right-hand diagram 
of Fig. 2. By checking the results with and without consid- 
ering the FF contribution, we find that there exist differences, 
where the obtained total cross sections of pp — > 7r"T with 
FF are suppressed compared with those without FF, which are 
shown in Fig. 2. In addition, with increasing £„«, the total 
cross section with FF goes down after reaching its maximum, 
while the total cross section without FF is goes up continu- 
ously. The forthcoming PANDA experiment can test these 
theoretical results. 
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FIG. 3: (color online). The total cross sections of pp — > n^rjc (left- 
hand) and pp — > n°J/if/ (right-hand) with different A values. Here, 
the black points with error bars are the E760 data [7]. 




FIG. 4: (color online). The center-of-mass frame angular distribution 
dcr/dCl of pp — > n'^Tjc and pp — » rfiJIij/. Here, the results are given 
by taking the range of £■„!! = 3.2-5.0 GeV or 3.4-5.0 GeV with step 
of 0.2 GeV for the rj^ or Jltj/ production. The diagrams in the first 
column are the results without FF [8] while the remaining diagrams 
are the results with FF, where we take A = 1 .9 GeV. The results at the 
lower and upper limits of £„„ range are highlighted with dot-dashed 
green and short-dashed red lines respectively, while long-dashed blue 
lines are the result at £„„ = 4 GeV. All results are normalized to the 
forward intensity. 



With pp — > nPrjc and pp — > jtPjI4i as example, we also 
present the variation of their total cross sections with different 
A values as shown in Fig. 3, where we take several typical 
values in the range of A = 1.3 ~ 2.3 with step of 0.2 GeV. 
These results show that the total cross sections of pp 
and pp — > n^jjxli depend on the value of A. 



n\. 



Besides providing the information of total cross section, we 
also give the result of the center-of-mass unpolarized angular 
distributions of pp — > Ti^r\c and pf> — > .tt'V/i^ just suggested in 
Ref. [8]. Considering the FF contribution, the corresponding 
angular distributions da-jd^ are shown in the second column 




B. The obtained result for the pp — > oj'V reactions 



FIG. 5: (color online). The center-of-mass frame angular distribution 
dcr/dO. of pp -^ ?rVcO, PP — » ^Vci ^nd pp — > n'^i//' corresponding 
to E„„ = 3.6-5.0 GeV, 3.8-5.0 GeV, 4.0 - 5.0 GeV with step of 0.2 
GeV, respectively. Here, the results without and with FF correspond 
to these diagrams in the first and the second rows, respectively. The 
style of the highlighted lines is arranged in the same way as in Fig. 
4. 



of Fig. 4, where the resuhs without FF are also listed. In 
addition, the center-of-mass unpolarized angular distributions 
of pp — » ttVcO, pp — > ■'I'Vfi ^nd pp -^ jf'if/' axe listed in Fig. 
5. We can find that the angular distributions of pp -^ tt"*? 
with and without FF are anisotropic just shown in Figs. 4 
and 5. Another common peculiarity of these results in Figs. 
4 and 5 is that these obtained distributions become forward- 
and backward-peaked with increasing £„«, which is consistent 
with the conclusion in Ref. [8]. 

The comparison of the results in Figs. 4 and 5 also indi- 
cates the different behaviors of the center-of-mass frame an- 
gular distributions dcr/dQ. with and without FF, especially for 
pp -^ 7T"i//', pp -> ttVcO and pfy -^ n^Xd- For pfy -^ ji^ric, 
its dcr/dQ. distribution with FF become more forward- and 
backward-peaked than that without FF when taking the same 
Ecm value. There obviously exists a node in its differential 
cross sections with and without FF when taking = 90° in 
the center-of-mass frame or f = u, which is totally different 
from the corresponding result without FF, where is defined 
as the angle between the light meson and proton. In fact, this 
behavior is also supported by analyzing the detailed analytic 
expressions of this process. For pp — > ifJI^, we notice that 
its differential cross section with FF is close to when taking 
Ecm - 5 GeV and 6 - 90° in the center-of-mass frame. The 
situations of pp — > tt^i^', pp — > ti^^XcO and pp — > Ji^^Xci are 
similar to that of pp — > n^Jjij/. These results further indicate 
that the FF contribution should be considered in studying the 
charmonium production plus a light meson at the low energy 
pp interaction. 

The above investigation of the center-of-mass unpolarized 
angular distributions of /9p — » .tt'^*!' can provide valuable infor- 
mation to the design of the PANDA detector and the analysis 
of the experimental data. What is more important is that fur- 
ther experiment can test the calculated total cross sections and 
the dcr/dQ. distributions of pp — > tt'^*!'. 
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FIG. 6: The total cross sections of pp — > wrjc (left column) and 
pp — » (jjJ/ifi (right column) dependent on £"„„. The diagrams in 
the first row are the result without FF, while the rest diagrams are 
the result with FF. Here, we take three different combinations of g^ 
and K^ to illustrate the total cross sections dependent on (g^, k^) just 
shown in (a)-(d). Among these four diagrams, diagrams (c) and (d) 
are obtained by taking A = 1.9 GeV. Diagrams (e) and (f) show 
the total cross sections with different A values and the typical value 
(g„,/.„)=(I2.2,-0.I2). 



In the following, we illustrate the results of pp -^ oj^. Just 
shown in Table. I, there exist different values of g^i and k^ for 
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FIG. 7: The predicted total cross sections of pp — > oj^ con'espond- 
ing to typical values A = 1.9 GeV and (g,„/^^) = (12.2,-0.12). 
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the ppo) coupling. With pp -^ (dtJc and pp -^ ajj/tfr as exam- 
ple, we list the variation of the total cross section of these two 
processes with Ecm when taking three typical combinations of 
go, and K^. Our calculation indicates that the total cross sec- 
tions of pp — > ajric and pp — > ojj/ip are dependent on the 
value of (gc^, Ka,) whether we consider the FF contribution or 
not (see Fig. 6 (a)-(d) for more details). The results shown 
in Fig. 6 (a)-(d) indeed indicate that the pp -^ ai^' processes 
can be applied to test the values of {go,, k^,) listed in Table. I. 
It is obvious that the total cross sections of pp — > corjc and 
pp — » coJIif/ are suppressed by FF. Since we cannot constrain 
the A value, in Fig. 6 (e)-(f) we discuss the dependence of the 
cross section of pp — > urjc and pp — * ajJIif/ on A. For other 
charmonium productions with o) meson, the behavior of their 
cross sections is similar to that of pp — > curie and pp — > aiJ/tfr. 
In Fig. 7, we further list the total cross sections for the 
pp — > a/i' processes taking typical values A = 1 .9 GeV and 
{gw,K:w) = (12.2, -0.12), where we take the same cutoff A as 
that of pp — > TTpJIip. Our calculation shows that the total cross 
section of pp — » utjc is the largest one among all charmonium 
productions discussed here. The production cross section of 
J I ill is roughly 50 times smaller than that of rjc, while the pro- 
duction cross section of XcO is 5 times larger than that of Jlif/. 
We also notice that the total cross sections of pp — > oj^ be- 
come stable after reaching up to their maximums with increas- 
ing Ecm, which is different from the situation of pp — > 7r°*i' 
discussed above. Just because of the considerable cross sec- 
tions of Tjc and XcO productions, pp — > utjc and pp — > uxm 
can be as the ideal channels to study the charmonium produc- 



FIG. 8: (color online). The obtained center-of-mass frame angular 
distribution dcr/dil of pp — > ojtjc and pp — > coj/if/. Here, the results 
are given by taking the range of fi'cm = 3.8-6.0 GeV or 4.0-6.0 GeV 
with step of 0.2 GeV for the ij^ or J/ifr production, respectively. The 
diagrams in the first and third rows are the results without FF while 
the remaining diagrams are the results with FF, where we take A = 
1.9 GeV. The diagrams with the same coupling constants (g^.^'w) are 
listed in the same column. The results at the lower and upper limits 
of £„„ range are highlighted with dotted green and short-dashed red 
lines respectively, while long-dashed blue lines are the result at £„„ = 
5 GeV. 



tion at PANDA. 

Additionally, with pp — > totjc and pp -^ coJ/tfr as example, 
we also give their center-of-mass frame angular distributions, 
which are shown in Fig. 8. There exists obvious difference 
among the results calculated with and without FF. We can do 
a cross check for our calculation, where we take A — > oo, 
which denotes that the ppio and pp^ interactions are treated 
as the point-like vertices. We find that the result with FF is 
gradually consistent with that without FF when A tends to 
infinity. We also find the node in the differential cross section 
with FF when taking 9 = 90° in the center-of-mass frame, 
which is similar to p/? -^ ;r"*I'. The investigation of ;9/? — » a/i' 
further shows that the FF contribution cannot be ignored in the 
calculation. 



C. Total cross sections for other processes 

We can easily extend the fonnulae of pp — » n^'V and pp -^ 
(jjV to study pp — > jfV and pp — > p^ respectively, where we 
only need to replace the corresponding coupling constants and 
masses. The coupling constant of the pprj interaction is taken 
as gpp^ - 11.5 [11]. In Ref. [22], the ppp coupling constants 
{gp,Kp) - (3.249,6.1) obtained by the Bonn full model. In 
addition, the masses of 77 and p mesons are M^ = 547.9 MeV 
and Mp = 775.5 GeV [3], respectively. 

In Fig. 9, the total cross sections of pp — > ifV and pp -^ 
pV are calculated by including the FF contribution, where we 
take the typical cutoff A = 1 .9 GeV. The result also shows 
that the total cross section of the ric production is the largest 
one among all charmonium productions by the pp — > rfV or 
pp — > p'V processes. 
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FIG. 9: The total cross sections of the processes pp — > rfV (left- 
hand) and pp — > pT (right-hand) dependent on £■„«■ Here, we take 
the typical value A = 1.9 GeV. 



IV. DISCUSSIONS AND CONCLUSION 

As one of the main physical aims of PANDA, studying the 
charmonium production at the low energy pp interaction is 
an interesting research topic not only for experimentalists but 
also for theorists. For investigating the production of char- 
monium plus a pion by the pp annihilation, the production 
mechanism shown in Fig. 1 was proposed in Ref. [5], where 
the pp — » TfJIijj cross section was calculated. According to 
this mechanism, the authors in Refs. [6, 8-10] further devel- 
oped the model of the production of charmonium plus a light 
meson and applied it to calculate the pp — » ;r"*I' processes, 
where the total cross sections and the unpolarized angular dis- 
tributions were obtained. Although there are many theoretical 
groups dedicated to the production of charmonium plus a pion 
by the pp interaction, further theoretical study of these pro- 
cesses is still valuable, where several further developments of 
the model were indicated in Ref. [10]. 

Considering the present research status of the production of 
charmonium plus a light meson by the low energy pp inter- 



action, in this work we revisit this interesting research topic. 
Different from the former work in Ref. [8], in our calcula- 
tion we consider the FF contribution to each interaction ver- 
tex, which was proposed as one of the possible developments 
of the model [8]. We calculate the total cross sections and the 
center-of-mass frame angular distributions of the pp — > if*^ 
processes. Since the E760 experiment measured the cross sec- 
tion of pp — > Tf'JIij/ [7], we compare our result with the E760 
data, which indicates that the calculated total cross section of 
pp — > tt"*!' overlaps with the experimental data, where the 
cutoff A in FF is taken as 1.9 GeV. The inconsistence between 
the experimental data and the result without FF in Ref. [8] 
are alleviated by considering the FF contribution, which in- 
dicates that the FF involved in each interaction vertex can- 
not be ignored. Adopting the same A value, in this work we 
also give other charmonium production cross sections and the 
corresponding center-of-mass frame angular distributions of 
pp — » tt"*!'. The difference of the results with and without FF 
further shows that we should consider the FF contribution to 
our calculation. Thus, these results obtained with FF can be 
served as further experimental investigation of charmonium 
production at PANDA. 

Besides studying the charmonium production plus a pion, 
in this work we also calculate the pp -^ u^ processes. The 
results of the total cross section and the center-of-mass frame 
angular distribution of pp -^ diV also indicate the distinct 
difference of the calculations with and without FF, which is 
similar to the situation of pp — > Tr"*!*. Different from the ppn 
coupling, the ppo) interaction is related to two independent 
coupling constants g^ and /<■„, which make pp — > cJ^ an ideal 
channel to test different theoretical values of g^, and k^,, where 
our results also show the total cross sections of pp — > cJV are 
dependent on the values of gca and /c„. 

In order to reflect the completeness of the study of the char- 
monium production plus a light meson, we also apply the for- 
mulae of pp — > tt"^ and pp — > (a^V to calculate the pp — > jfV 
and pp -^ p^, respectively. The predicted total cross sections 
of pp -^ rfV and pp — > p*? are accessible at the forthcoming 
PANDA. 

By the systematic investigation of the charmonium produc- 
tion plus a light meson at the low energy pp interaction, we 
notice a common behavior of the charamonium production, 
i.e., the ric production cross section is the largest one among 
all discussed processes. This fact shows that the low energy 
pp interaction is an ideal platform to produce //^ . In addition, 
we also find that the charmonium production cross sections 
satisfy the relation cr,^^ > cr^^ > crjj^ > ct"^^, > cr^', which 
does not depend on the associated light mesons. 

The authors in Ref. [6] calculated the total cross section 
of the //(/' and ifr' productions pp — > mJ/t//{iff') by relating 
these processes with J/iff{t//') — > mpp. Thus, we also make 
a comparison of the results obtained by us and those listed 
in Ref. [6]. As shown in Table. II, our result with FF of 
pp — > cokjj' ,p4i' are consistent with those given in Ref. [6] 
while the cross sections of pp — > 7:^i(/',rji(/' calculated by us 
are slightly smaller than those in Ref. [6]. The cross section 
of pp — > pJ/4r shown in this work also falls into the range 
of the predicted cross section of pp — > pJ/iff, where the up- 



TABLE II: The comparison of the total cross sections of pp — > 
mj/if/(tf/') obtained in this work (the second column) and given in 
Ref. [6] (the third column) and the corresponding £„„ value. Here, 
our results are the typical values of the cross section with FF when 
taking A = 1.9 GeV. 



Reaction 


a-pF (pb) 


0-™' (pb) [6] 


Ea„ (GeV) 


pp -^ n°JI^ 


116 


420 ± 40 


4.28 


pp -^ rjj/i// 


36 


1520 ± 140 


4.57 


pp — » U)J/tff 


156 


1900 ±400 


4.80 


pp -^ pJlip 


127 


<450 


4.80 


pp -^ ;r>' 


28 


55 ±8 


5.14 


pp -^ TJlj/' 


9 


33 ±8 


5.38 


pp —> dltf/' 


40 


46 ±22 


5.60 


pp -^ pi//' 


32 


38 ±17 


5.59 



per limit of the cross section of pp — > pJ/i^ was given in 
Ref. [6]. There also exist the differences of the results of 
pp —^ 7T^J/tl/,riJ/tl/,ojJ/tfr from this work and Ref. [6], i.e., 
total cross sections presented here are about 4 times, 42 times 
and 12 times smaller than those for pp — » t^'JI^', tiJI4'^ coJ/iff, 
respectively. 

In summary, in this work we systematically study the char- 



J 
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monium production plus a light meson in the pp interac- 
tion, and predict the corresponding total cross section and the 
center-of-mass frame unpolarized angular distribution, which 
provide valuable information to the experimental investigation 
of this kind of reaction. As an ideal experiment for studying 
the charmonium production, the forthcoming PANDA exper- 
iment can directly verify the prediction given in this work and 
test the charmonium production mechanism at the low energy 
pp interaction adopted here. 
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Appendix A: The difTerential cross section of pp — » (^iT 

By Eq. (5), we can easily deduce the expressions of the 
unpolarized differential cross sections, i.e., 






^OJ^T], 



PP^OJIJ, 



s(s — 4rnV)x-y^ 



[l-^((-/?2(r2 + 2)f- + if + IR" + 2r^R^ + 2R- f)xy) + K{3r'-R^f- - (if + 4r^R^)xy) 



do-\ 



dcrX 



--K^{{r^ + %)r^R"f - {r^f + Af - Ar^R'^ - 4rR~f + l6rR~)xy - 4(7?^ + f)x^y^)\, (Al) 

n . "fT, , , ((-(r^ + 2)(R^ - A)f + (2(r2 + 2)(r2 - 4) + Sr^ + 2R' - \2R^ + 16 + 2{R^ - 4)/ + f)xy) 

+l,K?f\{R^ - 4)/ + 2xy\ - -K^\?{? + 8)(/?- - A)f - {Ar^R^ + Ar^R" - l6r^R- + Ar\R^ - 8)/ 

8 

+(r^ - A)f)xy - Air" + R^ + /)xy]), (A2) 

^ "'-'^J/f'' h[ - (r2 + 2){R^ + 2)f + (2(r2 + 2)(R'- + 2) + 2(r^ + R^ + 2)/ + 2(r^ + R^f - 8 + f)xy 
s(s — Am~)x^y^ ' 

-2ji^/] + 2K[ir^(R^ + 2)f - (6r'* + lAp-R^ + 6r^f + f)xy] + -i?[ - r^{r^ + %){R^ + 2)/" + {2r^ + 6r'^R^ 

+ I6r'' + 2r*f + 32r^R- + 16r^f + r^f)xy + A{R- + f)x^f]], (A3) 



dcrX 



PP^OJXcl 



K " ^'' -, o (-7[ - R^ir' + 2)(R- - A)f + (2R^ + 6rR* - SR"^ - 16rR- + 2r'^R^ + 2R\r + R^ - 4)/ 

+{R^ + 2)f)xy - 2R^x^y^] - ^[ - 3r^R^{R^ - A)f + {2r'^R^ + 6r^R* - 2Ar^R^ - 2?R^f + {R^ + 2)f)xy 

R^ 

J- 
+2{r^ - 3/?2)xy] :^{r^R^{r^ + %){R^ - A)f + {2r^R^ + 2//?^ - 16//?- - 2ArR'^ + 96rR- + 2r^R^{r 

+A)f + (2r - r^R^ - 8)f)xy - 2{/ + r^R' + Ar^ + 2R'^ - UR^ + (r^ + 2R^ - 4)/)jcy - Ax^y^]], (A4) 



where a^ = g],p^l^n , a^v = glp^/^Jr , r = m^/m^ and R = 
m^/nip. As the dimensionless variables, x and y are defined 
as X = t/mjj - 1 and y = u/mi - 1, respectively. / denotes a 
dimensionless energy variable with definition f - (s - m^ - 

ml)/mj = -{x + y). 



Appendix B: The total cross section of pp —> wT 



The detailed formulae of the total cross section of pp 
w*}* are 



O-n 



' PP^i^Xai 



= ^-r^^^T-IT{[ - '^'('•' + 2)/j2 + (/ + 2R' + 2rR- + IR'fVx] + K[3r^R^fl2 - {if + ArR^)I,] 
s{s — 4mt,) * 

--/f-[(r- + d.)r^R^fl2 - {r^f + 4/- - Ar^R'^ - Ar^R-f + \6r^R-)Ix - \{R^ + /)Jo]}, 

-^\\ - (r^ + 2)(R^ - 4)fl2 + i2(r + 2)(r^ - 4) + Sr" + 2R^ - UR" + 16 + 2(R^ - 4)/ 



s{s — 4m-) ' 

+f)Ii] + 3Kr^f[{R^ - 4)/J2 + 2Ji] - -K^[r^{r^ + 8)(7;' - 4)fl2 - {AR'^R^ + Ar^R"^ - l6r^R^ 

8 

+4r\R^ - 8)/ + (r^ - 4)f)Ii - 4(r + R^ + /)Jo]}, 



O-pp^OjiJ/lf,"!) - 



Q'wQ'(y/iA">)'"p 
s(s — 4m?) 



[2[ - (? + 2)(R^ + 2)fl2 + {2{r^ + 2){R^ + 2) + 2(r^ + R^ + 2)/ + 2(r- + R^ - 8 



(Bl) 



(B2) 



1 0, 



o-„ 



+/-)Ji - 2 Jo] + 2/f[3r2(/?2 + 2)fl2 - (6/ + Ur^/?^ + 6r^f + f-)Ii] + -k-[- r{r + 8)(/?- + 2)/ J, 
+(2r* + 6//?2 + 16/ + 2rV + 32r27?' + 16rV + '"V^Ki + 4(7?^ + /)Jo]}, (B3) 

CK Ct ITl 

n " ^" '! {-r[ - /^^(r^ + 2)(/?^ - 4)/j2 + {2R'' + 6r^R^ - 8R'^ - 16r^R' + 2r^R^ + 2R\r^ + R^ - 4)/ 
s{s — 4mt) ^R-^ 

+{R- + 2)f)Ii - 2R^Io] - -^[- 3r^R^{R^ - 4)fl2 + {2r''R^ + er^R"^ - 24r^R^ - 2p-R^f + (R- + 2)f)Ii 

+2(r^ - 3R^)Io] 2^r^R^(r^ + 8){R^ - 4)fl2 + {2r''R^ + 2rV - 16/7?^ - 24r^R'^ + 96r^R^ + 2r^R^(p- 

+4)/ + {2r^ - r-R- - 8)f)Ii - 2(/ + r^R^ + 4p- + 2R'' - 12R^ + (r^ + 2R^ - 4)/)Jo - 4J_i]}, (B4) 

I 



dx {xyY 
I,n(xi) - I,„{xo) with y^-x-f. 



Xq 



and 



ml - 2EpE^ - 2ppPa 



Xl 



ml - 2EpEco + 2pppa 



J_i(x)= -(ix^ + lxA, 
Jo(x)= X, Ji(x)=iln(^^ 



Iiix) = -T In 



x + f\ 1/1 1 



P\x + f X 



When deducing these expressions of pp — > co^, the momenta 
and energies of proton and co meson are related to the the Man- 
delstam variables by 



1 / 2\l/2 S^'^ 



1 (m\, - mlY - 2(ml + m^)s + s^ 



^'" 5^172 (^ - '4 + '«^)- 
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